Quantum capacitance and density of states of graphene 



O 
(N 



CD 

S 

a 

o 
o 



o 

On 
O 

o 



X 



S. DroscherQ P. Roulleau, F. Molitor, P. Studerus, C. StampferQ and T. Ihn and K. Ensslin 
Solid State Physics Laboratory, ETH Zurich, 8093 Zurich, Switzerland 

We report on measurements of the quantum capacitance in graphene as a function of charge carrier 
density. A resonant LC-circuit giving high sensitivity to small capacitance changes is employed. 
The density of states, which is directly proportional to the quantum capacitance, is found to be 
significantly larger than zero at and around the charge neutrality point. This finding is interpreted 
to be a result of potential fluctuations with amplitudes of the order of f 00 meV in good agreement 
with scanning single-electron transistor measurements on bulk graphene and transport studies on 
nanoribbons. 



The density of states of a given quantum system is im- 
portant for understanding the electrical conductivity and 
other electronic response functions. Direct measurements 
of the density of states via the quantum capacitance have 
been successfully used, for example, on semiconductor 
samples containing a two-dimensional electron gas [l|, [1] 
and carbon nanotubes Q. Graphene sheets have been 
subject to theoretical studies of the quantum capacitance 
0, Q and first measurements have been carried out re- 
cently 01 . Techniques such as scanning tunneling mi- 
croscopy 3 were applied to probe the local density of 
states of graphene as well. We present experiments on 
a structure covered with a top gate similar to those in 
Ref. and interpret our results based on the presence 
of disorder in the graphene device. 

The samples were fabricated using mechanical exfo- 
liation of graphite powder and subsequent deposition 
onto 285 nm Si02 [9| on an undoped silicon substrate 
(see Fig. [Ija)). It was verified that the investigated 
layers consisted of a single atomic sheet using optical 
microscopy, AFM microscopy as well as Raman spec- 
troscopy [13, [iH ■ Ohmic contacts were patterned with 
electron beam lithography followed by metal evaporation 
of Cr/Au (2 nm/40 nm). In a second electron beam 
lithography step a mask for reactive ion etching was de- 
fined to structure the graphene sheet. The patterned top- 
gate consisted of Ti and Au evaporated onto an AI2O3 di- 
electric. The alumina dielectric was obtained in cycles of 
depositing 1 nm of aluminum followed by a 3 min period 
of increased oxygen pressure in the deposition chamber 
leading to the complete oxidation of the thin Al-film. Re- 
peating the cycle four times resulted in an oxide thickness 
of approximately 12 nm. A scheme of the cross section 
and an AFM image of the device investigated here are 
shown in Figs. [IJa) and (b). 

All measurements were carried out in a He-bath cryo- 
stat at a base temperature of 1.7 K. Besides transport 
measurements through the graphene sheet using stan- 
dard lock-in techniques, the capacitance between top gate 
and graphene sheet was recorded using an LC-circuit (see 
Fig. [He)). The oscillator was made of an inductor [L 
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FIG. 1: (a) Schematic cross section of a top gated graphene 
sample. The graphene flake (black) is contacted with gold 
electrodes (orange) and covered partly by an alumina/gold 
top gate (grey /red), (b) Atomic force microscope (AFM) im- 
age of the device studied here. The electrodes are colored 
corresponding to the scheme in (a) . (c) Circuit scheme of the 
measured system. 



= 100 /zH) placed at room temperature and the capaci- 
tance of the wiring inside the cryostat Ccabics ~ 340 pF 
in parallel to the sample capacitance Cg. The sample ca- 
pacitance can be decomposed into two components I, 2|: 
the geometric capacitance Cg with the quantum capaci- 
tance Cq in series. The circuit was driven with an oscilla- 
tion amplitude of 20 mV at its self-resonant frequency of 
850-900 kHz. Variations of the sample capacitance were 
measured as a change of this frequency with a relative 
accuracy of 2.5xl0~^, giving a sensitivity to capacitance 
changes of the order of 170 aF. 

Fig. [5] (a) shows a measurement of the resonance fre- 
quency /res as a function of applied top gate voltage. A 
maximum in frequency is observed at the Dirac point 
where Cq is lowest and therefore dominates Cg. As the 
density is increased Cq increases and its contribution 
to Cs becomes negligible beyond Cq = Cg. The con- 
stant background value of /res far away from the charge 
neutrality point is hence given by Cg -I- Ccabics- Since 
Ccabies IS large compared to Cg and the circuit indue- 
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FIG. 2: (a) Resonance frequency as a function of applied top 
gate voltage Vtg- Inset: Two-point resistance obtained in 
transport measurements, (b) Relation between Fermi energy 
and top gate voltage as deduced from eq. ((3]). The differ- 
ent curves (black lines) are obtained for varying values of Cg 
(from outermost to innermost curve: 4.8 fF //ira^ , 5.8 fF /fim^ , 
6.8 fF/^m^). The dashed curve shows the theoretically ex- 
pected dependence for a perfectly clean graphene sheet. To 
compensate for the shift of the Dirac point towards negative 
voltages the horizontal axis in (a) and (b) is offset by 0.3 V. 



tance L is known, Ccabics can be determined from the 
measured frequency far away from the Dirac point. Con- 
sidering the circuit model of our setup shown in Fig. [T] 
(c) we find that the change of the capacitance AC com- 
pared to the background is described by Cg/(Cg + Cq). 
From the resonance condition of such an LC-circuit AC 
can be determined by expanding /o for small variations 
leading to AC = -2(Ccabios + Cg)A///o. We define 
(5m = A/(Vtg)//o as the relative change in resonance 
frequency (see right scale in Fig. [21(a)). 

A two-terminal transport measurement is displayed in 
the inset of Fig. [5] (a). The resistance maximum is 
located at the same top gate voltage as the frequency 
maximum in the main figure, demonstrating that these 
independent measurements give a consistent charge neu- 
trality point. 

In order to transform the top gate voltage axis in Fig. [5] 
(a) into Fermi energy without using assumptions about 
the density of states, we consider the electrostatics de- 



scribing the structure. An applied top gate voltage in- 
troduces a difference in the electrochemical potentials of 
the top gate contact and the graphene ^q- The ef- 
fective potential difference is given by 



A*G - Mm = |e|VTG = 



const. 



(1) 



where Ep is the Fermi energy, Ug is the charge carrier 
density, d is the thickness of the dielectric layer and the 
constant term includes the work function difference of 
the two materials, which can be neglected for the analy- 
sis provided that it is gate voltage independent. The first 
term in eq. ([T]) is the chemical potential of the graphene 
sheet, whereas the second term describes the electrostatic 
potential difference between the two capacitor plates as 
obtained from solving Poisson's equation. The carrier 
density Ug is given by the density of states integrated 
from the Dirac point over all occupied energies up to Ep. 
Differentiating eq. ([I} with respect to E-p and identi- 
fying the geometric capacitance Cg/A — ee^/d and the 
quantum capacitance Cq(VTG) = yields 



TG 



dEp 



1 



gq(FTG) 

c. 



(2) 



Integrating the differential equation ([2]) over Vtg and 
including the relation between the capacitances involved 
and the measured frequency leads to an equation for the 
Fermi energy 



lei f"^"^^ 

Sf(^tg) = — / <5,„(KfG)dKfG 



(3) 



where Vb is the gate voltage at the Dirac point (here 
V) and a = Cg/2Ccabios- 

As the geometric capacitance can not be deduced from 
the data it is left as a free parameter that has to be esti- 
mated by other methods. The plate capacitor geometry 
of the device gives Cg « 6 fF / ixvo? . Hall measurements on 
a bilayer Hall bar structure located on the same chip were 
used to deduce the charge carrier density and from this 
the geometric capacitance. A value of Cg k, 5.8 fF //xm^ 
is obtained by this additional method and hence an esti- 
mate is found for the demanded parameter. 

Fig. [5] (b) displays the dependence of E-p on Vtg for 
three different values of Cg. Assuming a perfectly clean 
graphene sheet, one would obtain a relation indicated by 
the dashed curve. Close to Vtg = a linear behavior 
is observed leading to a large increase of Ep for small 
changes of Vtg which levels off towards the boarders 
of the investigated bias window indicated by the dot- 
ted vertical lines. The error bars show the uncertainty of 
the values calculated from eq. ([3]) assuming a Gaussian 
probability density distribution for (Sm with the width 
(7=2.5 xlO~^ determined from the variance of the noise 
seen in Fig. [2] (a). With increasing Ep the error grows 

with Ep (Vtg ~ ) • 

With the result in Fig. [2] (b) at hand, allowing us to 
convert the Vtg axis in Fig. [11(a) to Ep, we now proceed 
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FIG. 3: Density of states as a function of Fermi energy. The 
black solid lines show the experimental data assuming dif- 
ferent Cg (from outermost to innermost curve: 4.8 fF//im^, 
5.8 fF//im^, 6.8 fF//jm^). Error bars are indicated and the 
theoretically expected density of states for a perfectly clean 
graphene sheet is drawn as the dashed line. The bias window 
defined in Fig. [2] is plotted with dotted lines. 



with the aim of transforming the A/ / /o axis to a density 
of states. The quantum capacitance and therefore the 
density of states is related to the geometric capacitance 
and the change of the total capacitance by 

C,^C,^-^. (4) 

Om 

However, the measurement uncertainty gets very large 
outside the interval indicated by the vertical dashed lines 



in Fig. [2] because the signal to noise ratio quickly ap- 
proaches zero. A careful analysis of the data using the 
methods of Ref. 12 leads to the density of states V{Ep) 
plotted in Fig. [3] for different Cg between 4.8 fF / 
and 6.8 fF / ixir? . This is the main result of this paper. It 
shows a flat bottom in the energy interval Ey = ±50 meV 
at about 1x10^^ eV~^m^ with the tendency to increase 
outside this interval. Error bars indicating the FWHM of 
the underlying probability density function for TD^E-p) are 
shown for selected data points. The dashed line in Fig. [3] 
shows the theoretical density of states of ideal graphene 
given by T>{Ey) = 2Ey /t^{vyK)^ for comparison. [l3j 

The data demonstrates that, in contrast to theory, 
there is a finite number of states close to the charge 
neutrality point. We attribute this observation to the 
presence of local potential variations within the graphene 
sheet on the energy scale of about lOOmeV. From trans- 
port measurements on nanoribbons p^ - [l8j and studies 
using a scanning single-electron transistor [§1 a compa- 
rable magnitude of disorder was inferred. This suggests 
that the gate oxide does not induce additional disorder 
on this energy scale. 

In conclusion, we have demonstrated a highly sensitive 
method for capacitance measurements using a resonant 
circuit. We have extracted the density of states of a single 
layer graphene sheet in which the density was tuned by a 
top gate. The magnitude of the disorder potential of the 
device was determined and is in good agreement with re- 
sults obtained in other studies. However, in order to map 
features like Landau levels or the one-dimensional den- 
sity of states of graphene nanoribbons [El the potential 
fluctuations have to be reduced drastically by improving 
the sample quality. 

We thank K. von Klitzing, J. Martin, M. Gcllcr and 
I.V. Zozoulenko for helpful discussions. 
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